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ABSTRACT 

Space-based gravitational wave interferometers are sensitive to the galactic population 
of ultra-compact binaries. An important subset of the ultra-compact binary popula- 
tion are those stars that can be individually resolved by both gravitational wave in- 
terferometers and electromagnetic telescopes. The aim of this paper is to quantify the 
multi-messenger potential of space-based interferometers with arm-lengths between 1 
and 5 Gm. The Fisher Information Matrix is used to estimate the number of binaries 
from a model of the Milky Way which are localized on the sky by the gravitational 
wave detector to within 1 and 10 deg 2 and bright enough to be detected by a magni- 
tude limited survey. We find, depending on the choice of GW detector characteristics, 
limiting magnitude, and observing strategy, that up to several hundred gravitational 
wave sources could be detected in electromagnetic follow-up observations. 
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1 INTRODUCTION 

A variety of detector concepts for space-based gravitational 
wave interferometers have been proposed, th e most well 
studied concept being LISA (|Bender et"alll998l ). It was un- 
derstood early on that the most numerous source class radi- 
ating in the band covered by LISA-like detectors will be the 
galactic population of ultra-compact binaries (UCBs) com- 
prised of pairs of stellar remnants: white dwarfs, neutron 
stars or black holes. The gravitational radiation from these 
UCBs will be the dominant signal in the frequency band 
covered by LISA-like detectors. 

Early estimates of the composit e signal from the UCB s 
|Evans et all Il987l : iHills et all Il99d ; iHills fc Bender] 119971 ) 
demonstrated that the signals of the vast majority of the 
galactic binaries will overlap and be unresolvable from 
one another, forming a limiting foreground (or "confusion 
noise") for space-based gravitational wave detectors. Later 
studies based on popula t ion synthesis ([Nelemans et a] 200 ll: 
iBenacquista et a 1 2004; E dlund et aj 120051 : iTimpano et all 
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l2006l : iRuiter et all |2006| ) have borne this expectation out. 
Detailed data analysis studies have shown that ~ 10 4 
individual binaries could be resolved out of the fore- 
ground by a gravitational wave observatory like LISA 



ground by a gravitational wave observa tory like blbA 
{Timpano et all l2006l: ICrowder fc Cornishl [20071 : iLittenberd 
l201ll : lNissanke et alH2012l ). 



A subset of the resolvable binaries will be detectable 
electromagnetically. The purpose of this work is to assess 
the multi-messenger potential for different space-based de- 
tectors spanning the trade-sp ace of future mission designs. 
This builds off previous work (|Coorav et al| [2003: Nel emans] 
120061 , 12009] ) demonstrating the feasibility of follow-up obser- 
vations for high-frequency UCB sources. We estimate the 
total number of multi- messenger sources by beginning with 
a pop ulation synthesis model of the galaxy {Nelemans et all 
|2004| ). complete with optical magnitudes. From this we pro- 
duce a magnitude limited source catalog, then estimate how 
well each system will be localized on the sky by different 
gravitational wave detector configurations. Using hundreds 
of Monte Carlo realizations over the spatial distribution of 
the galaxy and the UCB orientations, we find tens to hun- 
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Table 1. Gravitational wave detector configurations used in this 
study. Configuration 1 corresponds to eLISA. Configuration 5 is 
the classic LISA design. All simulations were for two year mission 
lifetimes. 



dreds of sources that can be observed both electromagneti- 
cally and gravitationally. 

The information encoded about the UCBs in each of 
the two spectrums is highly complementary, enabling tests 
of general relativity, full measurement of the physical pa- 
rameters enabling constraints on binary synthesis channels, 
and new methods of probing the close interaction dynamic s 
of the compact stars i|Cutler et a]||2003l ; Istroeer et alll2005l ). 



2 DETECTORS 

For a gravitational wave observatory, the limiting sensitivity 
as a function of frequency is dominated at low frequencies by 
acceleration noise S a , while the "floor" , where the detector is 
most sensitive, is dominated by position measurement noise 
S x . Table [1] contains the parameters used for the detector 
configurations in this study. These parameters can be used 
to compute the noise power spectral density 

Sn(f) = S ga i(/) + (4/3)sin 2 M[(2 + cosu)& 

+ 2(3 + 2cosu + cos2u)S' [l /(27r/) 4 ] , (1) 

where u = 2nf£/c and S ga i is the contribution to 
the ins trument noise from the unresolved UCB fore- 
ground (|Timpano et alll200r3 ). 

The configurations we will highlight correspond to the 
classic LISA design (£ — 5 Gm), as well as two shorter arm- 
length configurations (1 = 1 Gm and £ = 1 Gm) in order 
to cover a variety of plausible mission configurations. The 1 
Gm configuration is similar to the eLIS A mission being con- 
sider ed by the European Space Agency (|Amaro-Seoane et all 
2012). We use an observation time of two years for each con- 
figuration. 

This suite of detectors provides a broad palette to il- 
lustrate the observational capabilities of these instruments 
with regards to the UCBs. A classic depiction of the per- 
formance for these interferometers is a plot of the average 
sensitivity curve in strain spectral density versus frequency 
(ILarson et a j |200d ). as shown in Fig. [1] The eLISA con- 
cept is the only one which uses a 4-link configuration. The 
doppler ranging between each spacecraft in the constella- 
tion is accomplished using two laser links. Thus the 4-link 
design is a single- vertex interferometer, while the 6-link de- 
signs allow for three (coupled) interferometers. This differ- 
ence accounts for an additional improvement in the 6-link 
sensitivity curves by a factor of ~ V2 at frequencies where 
the UCBs are found. 
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Figure 1. Sensitivity curve for each of the detector configura- 
tions in Table [T] The solid lines show the sensitivity set by the 
measurment noise while the dashed curves include an estimate 
of the UCB confusion-limited foreground. Over-plotted are the 
brightest UCBs in our simulated catalog (green crosses), and the 
known verification binaries (blue stars) 

3 DISCOVERING NEW VERIFICATION 
BINARIES 

The focus of this work is to study the population of de- 
tectable UCBs in the context of multi-messenger astronomy. 
We will focus on the sources detected via GWs which could 
potentially be identified electromagnetically. There is a sep- 
arate class of UCBs, the "verification binaries," which are 
known low-frequency GW sources with AM CVn serving 
as the archetype. There are ~ 30 known verification bi- 
naries, ~ 5 — 10 of which could be identified by the GW 
detec tors considered here, with sources still being discov - 
ered (|Nelemandl2"oi"ll ; iRoelofs et al2007l ; lBrown et aj|201ll ). 
This study does not include the known verification binaries 
in the galaxy catalogs. Furthermore, many of the AM CVn 
systems would not be localized well enough by the GW mea- 
surement alone to warrant simple electromagnetic follow-up 
observations. 

3.1 Binary selection 

The UC B population model i s essentially identical to that 
found in lNelemans et all (|2004l ). so the positions and ages of 
the systems are based on the Boissier & Prantzos (1999) 
Galactic model. We use the white dwarf cooling tracks 
based on Han s en (1 999) as shown in the Appendix of 
iNelemans et~al (|2004t ). We convert the luminosities to V- 
band magnitudes using zero-temperature white dwarf radii 
and simple bolometric corrections based on the effective 
temperature. This should suffice for this initial estimate of 
the potentially detectable population, but can be improved 
using detailed WD cooling models in the future. We deter- 
mine the absorption as in INelemans et 3 |2004h based on 
the Sandage (1972) model, but correcting for the fact that 
the dust is more concentrated than the stars, so we use 120pc 
as scale height for the absorption. 

To construct the magnitude limited catalog, we begin 
with the entire binary population in the synthesized galaxy. 
The limiting apparent magnitude of a telescope is a function 
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of the aperture D, the exposure time t, and the properties 
of the detector use d for imaging and photometry |Schaeffer] 
1 199d ; [Howelll Il98^ ) . A rudimentary fit to the limiting mag- 
nitude m using a telescope of aperture D (in m) and for ex- 
posure time t (in seconds) is given by m = 19.6 D 073 t 025 . 
Using commercial CCD detectors, a D — 0.5 m telescope will 
reach a photometric magnitude m — 21 in t ~ 75 s, where 
as a D = 1.0 m telescope will reach the same magnitude 
in t ~ 20 s. This paper examines the role of small to large 
aperture telescopes by examining a broad range of limiting 
magnitudes; lower bounds of m = 18 — 24 were chosen as the 
electromagnetic cutoff. All sky survey instruments such as 
LSST could further improve the number of candidates. The 
single exposure limit for LSST is expected to be m ~ 24, 
whereas the magnitude limit o f the final stacke d image is 
expected to be around m ~ 27 jlvezik et a]||201ll '). 



3.2 Gravitational wave detector response 

From the magnitude-limited catalog, we determine the num- 
ber of "bright" UCBs that will be well measured by the 
GW detector. To do so we must first estimate the confu- 
sion noise for each configuration. The instrument response 
to the galactic foreground is constructed by generating 
and co-adding waveforms for each source in the full sim- 
ul ated galaxy catalogue using the fast-slow decomposition 
in lCornish fc Littenbera ((2003). The confusion noise, S ga i, 
is empirically determined from the simulated data by it- 
eratively removing sources brighter than a running esti- 
m ate of the backgroun d. This procedure is first discussed 
in iTimpano ei~ail 120061') with an im proved implementation 
used here as in iNissanke et~al (|2012T ). 

With the confusion noise incorporated into the detec- 
tor sensitivity curves, we determine how well the GW de- 
tector can measure the source parameters of a UCB wave- 
form, using the we l l-worn Fisher Information Matrix I\j 
(| Cutler fc Flan agan 1994), the inverse of which approxi- 
mates the covariance matrix. There is no shortage of lit- 
erature highlighting short-comin gs of the Fisher to approxi- 
mate GW parameter errors e.g., Vallisneri (2008). However, 
given the scope of the problem we are addressing (hundreds 
of Monte Carlo's of thousands of detectable binaries) more 
rigorous parameter estima tion studies would be impractical 
(recently I Vallisneri (|201ll ) has proposed a way around this 
dilemma). On the other hand, the UCBs in which we are 
most interested - those that can be well localized on the 
sky - have atypically high signal to noise ratio, where the 
Fisher provides a good esti mate of the true parameter errors 
(ICrowder fc Cornish|[2007l '). 

For a binary to be considered "well localized" we require 
that the 63% confidence interval of the sky-location poste- 
rior distribution function subtends an area on the celestial 
sphere below some threshold dQ,. To bracket the capabilities 
of ground-based optical telescopes, we perform the analysis 
with dft ^ 1 and ^ 10 deg 2 . We estimate the area of the sky- 
location error ell ipse using the full cova riance matrix found 
by inverting r 4j (|Lang fc Hughesll2008l 'l. 

The number of well localized, bright binaries is com- 
puted for hundreds of realizations where we Monte Carlo 
over the orientation of each binary, as well their location 
within the Galaxy. For the orientation, we draw the inclina- 
tion i from a uniform distribution cos i = U[— 1, 1], and the 



polarization angle ijj and initial phase <p from U[0,27r]. We 
find up to a several hundred GW sources will be viable can- 
didates for electromagnetic follow-up searches, depending on 
the depth of EM survey and the GW detector characteristics 
(See the left-hand panel of Figs. [2] and [3]). 

3.3 EM detection strategies 

We now consider how to select candidates for follow-up ob- 
servations from the full GW catalog. Pointing telescopes at 
all of the GW sources localized within the adopted thresh- 
old would be inefficient, as we find between 10 3 and 10 4 GW 
sources in the full catalog will meet the dfl ^ 10 deg 2 thresh- 
old, while < 10% are likely to be brighter than m — 24, and 
only < 1% pass the m ^ 20 cut. 

Additional considerations need to be made to increase 
the efficiency of follow-up observing campaigns. We illus- 
trate two simple ways to isolate the GW sources that may be 
electromagnetically observable. These suggestions are sup- 
ported by calculations shown in Table [2] 

First, the large majority of UCB sources are confined 
within the galactic plane. Conversely, the magnitude lim- 
ited catalogs sample the local galaxy, which is much more 
uniformly distributed on the celestial sphere. Therefore, as 
a rough cut on the GW catalog, any binaries that are well 
localized but out of the galactic plane are good candidates. 
These are additionally attractive sources, as there will be 
less optical background and extinction against which the 
observing campaign will have to compete. We find between 
~ 20% and ~ 50% of the well-localized binaries in the 20 th 
to 24 th magnitude-limited catalogs have galactic latitudes 
\b\ 20°, while that fraction is reduced to ~ 1% for the full 
GW catalog. A uniform distribution of stars on the celestial 
sphere would have 66% of the stars with |6| ^ 20° 

The other strategy for identifying optical counterparts 
relies on estimates of the distance to the galactic binary. 
Typical UCB sources will undergo very little evolution of 
their orbital period during a space-borne GW detector's life- 
time. Without measurement of the rate of change of the 
gravitational wave frequency / the GW observation only 
constrains the overall amplitude of the signal without decou- 
pling the chirp- mass and the l uminosity distance d^ (|Schutj 
ll986l ; IStroeer fc VecchicboOrj l. For ~ 10-20% of the multi- 
messenger sources we sufficiently constrain / and A to mea- 
sure dL to within 20%, but astrophysical effects such as tides 
may impact the orbital evolution and thus bias the distance 
estimate. For the remaining systems in the GW catalog, we 
can use reasonable priors on the mass and mass ratio of 
white dwarf binaries to put meaningful constraints on dL 
from the amplitude measurement alone. 

Using only the amplitude, frequency and priors on the 
masses constructed from the population synthesis simula- 
tion, we find that the distribution of the most likely (ML) 
luminosity distances dl/ML is strongly peaked between and 
8 kpc - the distance to the galactic center - for the mag- 
nitude limited catalogs. The dLiviL distributions for the full 
well-localized catalog with no magnitude cut is more uni- 
formly distributed over a larger range. 

Our final consideration pertains to the expected optical 
light curves for UCB systems identified in the GW catalog. 
The population synthesis galaxy in our study is restricted 
to detached white dwarf binaries, as opposed to interacting 
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Figure 2. Number of binaries with sky- location resolved to within 1 deg 2 for each configuration as a function of limiting magnitudes. 
The left panel shows the total number of candidates, the right panel shows the subset of eclipsing binaries. The error bars represent the 
full range after Monte Carlo'ing over the location and orientation of each UCB system. Even the modest detection abilities (magnitudes 
m ~ 19) of small aperture telescopes can yield several electromagnetic counterparts; larger telescopes with deeper magnitude grasp will 
have significantly more sources that can be surveyed. 
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Figure 3. Same as Figure [21 except here we use an angular resolution threshold for the GW detector of dQ ^ 10 deg 



AM CVn systems. Without mass transferring from one star 
to the other in the binary, photometric variability is not 
guaranteed. The systems in the GW catalog that are best 
constrained are typically those at the high-frequency end of 
the population. This is to our advantage, because the shorter 
period binaries have a higher probability of eclipsing one 
another during an orbital cycle. 

We can put an additional cut on our EM/GW catalog 
by requiring the binaries to be eclipsing. (See the right-hand 
panel of Figs. [2] and El). From simple geometrical arguments 
|Coorav et al|[2003r ) the minimum inclination angle with re- 
spect to our line of sight that will produce eclipsing light 
curves is 

cos(t min ) ~ 0.3(//3.5 mHz) 2/3 (2) 

assuming all binaries have mass Mtotal ~ 0.5 M@ and radius 
-Rwd ~ 10 4 km. If we only consider binaries in the multi- 
messenger catalog with inclination angle less than t m i n , we 
reduce the total number of candidates by a factor of ~ 3. 



Nevertheless, we still find upwards of ~ 100 candidates for 
the large GW detector configurations and deep, wide field, 
optical surveys. Requiring eclipsing light curves significantly 
degrades the multimessenger potential for the 1 Gm configu- 
ration using catalogs limited to 20 th magnitude and dimmer 
- such EM follow-up surveys could come up empty. 



4 DISCUSSION 

We conclude that space-based gravitational wave detectors 
will be useful observatories for discovering new UCBs in the 
galaxy that could be observed electromagnetically, though 
deep, wide field, optical surveys may be required to produce 
large catalogs. We reach this verdict by considering a range 
of plausible near-future space-based gravitational wave de- 
tector concepts, and assess their measurement capabilities 
for magnitude limited catalogs of UCBs. Magnitudes for the 
constituents of each binary were derived from the population 
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|6| > 20° 


< 20% 
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1700 (370) 


0.02 (0.03) 
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0.08 (0.17) 


444 (197) 


0.22 (0.22) 
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0.14 (0.26) 
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signals. Here 


we enumerate 


the fraction 



of binaries that will make good candidates for electromagnetic follow-up observations. Plain numbers in the table correspond to the 
d£l ^ 10 deg 2 , while those in parenthesis correspond to the 1 dcg 2 threshold. Wc tabulate the average number of binaries N that meet 
the sky-resolution requirements (column 2), and then the fraction which are significantly above the disk (|6| > 20° - column 3), or have 
d,L measured to within 20% (column 4), the idea being that near-by binaries are more likely to be optically detectable - proximity can 
be inferred by |b| and/or determined through d^. The columns then repeat for the m ^ 20 and m ^ 24 magnitude-limited catalogs. 



synthesis simulations, and the gravitational wave measure- 
ment capabilities were estimated using the Fisher Informa- 
tion Matrix. Any UCBs that were brighter than our cho- 
sen magnitude limits (18-24) and located on the sky by the 
gravitational wave detector to within angular resolution dQ 
were considered multi-messenger candidates. We estimated 
the multi-messenger catalog sizes for both dfi ^ 1 and 10 
deg 2 . 

At the pessimistic end, we consider magnitude 18 lim- 
ited catalogs, and single-vertex interferometers with 1 Gm 
arm-lengths. The best scenario considered the classic LISA 
design and an optical telescope limited at 24* magnitude. 
The number of multi-messenger candidates was anywhere 
from a few to several hundreds over that range of detec- 
tor capabilities. If we put on the additional constraint that 
the sources must be eclipsing to allow for electromagnetic 
observation the counts were reduced by a factor of ~ 3. 

While most of the known verification binaries are AM 
CVn type stars, our study only considered detached white 
dwarf binaries, thus providing a very complimentary catalog 
of UCB multi-messenger systems. 

This work considered a conservative approach to find- 
ing multi-messenger UCBs, with competing criteria that 
strongly affect the expected population of systems de- 
tectable in both spectrums. Electromagnetic detections are 
most strongly affected by the magnitude limit of the detec- 
tion survey, a function of telescope aperture. By contrast, 
the gravitational wave detection catalogs of UCBs are ex- 
pected to have thousands of systems in them; most will be 
too faint to be detectable by any electromagnetic survey. 
However the gravitational wave localization criterion is a 
strong constraint on the multi-messenger catalog. We find 
that wide-field surveys (cft7 ^ 10 deg 2 ) yield more candi- 
dates than more narrow fields of view {dQ, < 1 deg 2 ) by 
50-100% for the full catalogs, and by a factor of 2-4 for the 
eclipsing binaries. 

We have estimated the number of UCB multi-messenger 
candidates without considering what could be done with 
joint GW and EM observations. Our follow-on effort will 
consider the science yield from joint observations of both 
the known verification binaries - mostly mass-transferring 
systems - and the close, detached binaries that will be dis- 
covered by space-borne gravitational wave detectors. 
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